Dynamics of holographic vacuum energy in the DGP model 



o 
o 

o 

CD 

Q 



Xing Wu/ Rong-Gen Cai,^ and Zong-Hong Zhu^'B 

'Department of Astronomy, Beijing Normal University, Beijing 100875, People's Republic of China 
"^Institute of Theoretical Physics, Chinese Academy of Sciences, 
P.O. Box 2735, Beijing 100080, People's Republic of China 

We consider tlie evolution of the vacuum energy in the DGP model according to the holographic 
principle under the assumption that the relation linking the IR and UV cut-offs still holds in this 
scenario. The model is studied when the IR cut-off is chosen to be the Hubble scale , the 
particle horizon _Rph and the future event horizon Rc\i, respectively. And the two branches of the 
DGP model are also taken into account. Through numerical analysis, we find that in the cases of 
in the (-f) branch and i?eh in both branches, the vacuum energy can play the role of dark 
energy. Moreover, when considering the combination of the vacuum energy and the 5D gravity 
effect in both branches, the equation of state of the effective dark energy may cross —1, which may 
lead to the Big Rip singularity. Besides, we constrain the model with the Type la supernovae and 
baryon oscillation data and find that our model is consistent with current data within la, and that 
the observations prefer either a pure holographic dark energy or a pure DGP model. 
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I. INTRODUCTION 

Recent SNe la and WMAP observations P, [3] have 
indicated that our universe is currently undergoing an 
accelerating expansion, which confront the fundamen- 
tal theories with great challenges and also make the re- 
searches on this problem a major endeavor in modern 
astrophysics and cosmology. The origin of the cosmic ac- 
celeration is still a mystery and is referred to as the dark 
energy problem. Various models have been proposed to 
solve this problem. They generally fall into the following 
two ways. One is to add an exotic energy component 
with negative pressure, that is, the dark energy, to the 
total energy budget of the universe. Among others, the 
most competitive candidate of dark energy sofar is the 
cosmological constant due to both its theoretical simple- 
ness and its great success in fitting with observational 
data, although it suffers from the cosmological constant 
problem 3 . Such problem is expected to be solved or al- 
leviated in the models of dynamic dark energy (see [3] for 
a more detailed review), which generally contains a scalar 
field evolving in time and driving the acceleration, just 
like the scalar field introduced for the inflation stage at 
early universe. In fact, the cosmological constant prob- 
lem is essentially a problem of quantum gravity. In quan- 
tum field theory, where the effect of gravity is neglected, 
the vacuum energy is determined by the UV cut-off kc, 
that is PA cx A:^. No matter how we choose the UV cut- 
off, be it the Planck scale lO^^GeV or the electroweak 
scale TeV, the value predicted by theory is far greater 
that that observed IQ-'^^GeV. bmce we are concerning 
problems at the cosmological scale, however, we should 
have take into account the effect of gravity. It is expected 
that the value of the cosmological constant or the vacuum 



energy would be predicted correctly from a complete the- 
ory of quantum gravity, which is still being explored. But 
at present the holographic principle, which is believed to 
be an important feature of quantum gravity, may shed 
some light on solving this problem. Follow the line of 
the holographic principle, the holographic dark energy 
model[7] is a promising candidate for solving the dark 
energy problem. In this model, the vacuum energy is 
no longer a time- independent constant, but evolves with 
time according to the holographic principle. The vacuum 
energy is related to the length measure on cosmological 
scale 



PA 



(1) 



where pA is directly related to the UV cut-off, L is the 
IR cut-off, Mp denotes the Planck mass and c is a nu- 
merical factor by convention, which is the parameter of 
the model. This relation of the entanglement of UV/IR 
was first proposed in fs'] , where L was first chosen as the 
Hubble scale. Then Hsu 9] pointed out that this would 
not lead to the desired equation of state. Finally, LiQ 
proposed the holographic dark energy model where L is 
the event horizon. And this model fits very well with 
current observations [l6j. 

The other way to solve the dark energy problem is 
to modify the theory of gravity at large scale, without 
resorting to any new energy component. For example, 
the f{R) theoryQ modifies the standard Einstein-Hilbert 
action to introduce an effective dark energy component 
in the Einstein frame. Here we focus on the DGP model 
Q , which describes our universe as a AD brane embedded 
in a 5D Minkowski bulk and explains the origin of the 
dark energy as the gravity on the brane leaking into the 
bulk at large scale. The model is described by the action 
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where gab is the bulk metric and /i^i, is the induced metric 
on the brane. The first term contains the 51? Ricci scalar 
whereas the second term contains the AD Ricci scalar 
on the brane, which is an extra term due to quantum 
effects, in contrast to the Randall-Sundrum scenario 11 1. 
The third term represents matter localized on the brane. 
And Sgh is the Gibbsons-Hawking boundary term. 

In this paper, we assume that the relation Eq.([T|) 
still holds in the DGP model, and consider the evo- 
lution of the vacuum energy on the brane (or the 
brane tension) according to the holographic princi- 
ple. Note that there are other models also gener- 
alizing the standard DGP model by adding a cos- 
mological constant (LDGP) [13], a Quiessence perfect 
fluid(QDGP)[ll, a scalar field(SDGP)[i3, or the Chap- 
lygin gas(CDGP) [iBl . Although the holographic dark en- 
ergy model is well consistent with observational data, it 
should be noted that the core of the holographic principle 
is that it relates the UV and IR cut-offs of a local quan- 
tum field system, which reflects some feature of quantum 
gravity. Thus even if the holographic vacuum energy had 
not played the role of dark energy, it would still be of sig- 
nificance to study on this problem. 

The paper is organized as follows. In section [Til we 
present the model under three cases of choosing the IR 
cut-off for the vacuum energy as respectively the Hubble 
scale -ff^^, the particle horizon i?ph and the future event 
horizon i?oh- We study the equation of state (EoS) of the 
vacuum energy and the effective EoS due to the combined 
effect of both the holographic vacuum energy and the 5D 
gravity effect. In section Hill wc use recent observational 
data to constrain the model and fit its parameters. We 
conclude this paper in the final section. 



II. THE MODEL 



which the branch is the self-accelerating solution in 
which the universe may enter into an accelerating phase 
in late time by virtue of pure 5D effect of gravity, while 
the (— ) branch cannot undergo an acceleration without 
additional dark energy component. Here the vacuum en- 
ergy is added in p = pm + Pa- Obviously the dynamics 
is typically different from the standard FRW cosmology. 
In addition, we have the usual equation of conservation 



P + 3H{1 



0, 



(5) 



where w is the equation of state (EoS). Here we assume 
that there is no interaction between matter and vacuum 
energy. Therefore both components obey the equation of 
conservation respectively and, in particular, for matter 
we have p,„ — pmoa^^, the same as that in usual FRW 
cosmology. We note that here the EoS of the vacuum 
energy evolves with time due to the holographic principle 
as shown later, as opposed to ACDM where wa = — 1. 

When applying the holographic principle in cosmology, 
a crucial problem is how to choose the IR cut-off L. In 
usual FRW sosmology, it is shownQ that only in the case 
of choosing L as the future event horizon can vacuum en- 
ergy play the role of dark energy. In the DGP framework, 
we consider L as the Hubble scale H^^ , the particle hori- 
zon i?ph and the future event horizon i?oh respectively. 
And then we investigate the evolution of the EoS of the 
vacuum energy. 

Intuitionally, we may guess that for one thing, in the (-) 
branch where there is no self-acceleration, adding a com- 
ponent whose EoS > — 1/3 will not lead to an accelera- 
tion whereas adding a dark energy component may cause 
the universe in this branch to accelerate; for the other 
thing, in the (-I-) branch, no matter whether the holo- 
graphic vacuum energy we add can in itself play the role 
of dark energy, the combined effect may lead to an accel- 
eration due to the contribution from the self-acceleration 
of this branch. 



We assume a fiat, homogeneous and isotropic brane in 
accordance with the result of the WMAP observation^. 
Following ^ , the Friedmann equation is 



or equivalently 



3M2 
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where H = a/ a is the Hubble parameter and rc = 
Mp/2M|g^ is the distance scale refiecting the competi- 
tion between 4D and 5D effects of gravity. For ^ Vc 
(early times), the 4D general relativity is recovered; for 
> Tc (late times), the 5D effect begins to be sig- 
nificant, e = ± represents two branches of the model of 



A. L as H-^ 

By Eq.lT]), the vacuum energy is p\ = Sc^MpH^. 
For convenience we insert this result into the Friedmann 
Eq.gl). After defining f7„ = g^^ff^ = f7„o(l + zf, 
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where Hq is the Hubble parameter at redshift z = 0, the 
above equation can be transformed into 



(1 - c')i^r ^ 2e^^^ 



n^a = 0, (6) 
Then we can solve the above equation to get E ~ H/Hq. 



Case 1: c = 1 
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where e = +1 corresponds to the contracting sohi- 
tion, and e = — 1 represents the expandmg solution 
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the conservation equation of oc a ^ there- 

fore Pa oc oc (X a^^, namely, the vacuum 



In this case, the solution is easy to find by virtue of 

m ■ Pm ^ ^ 

energy decreases faster than pm and it cannot be 
dominant in late time. Thus we do not consider 
this case as of physical interest. 

When c 7^ 1 we write the general solution 



1-C2 



(9) 



• Case 2: c > 1 Since flr^ is constant, the part under 
the square root on the RHS of Eq. ^ may become 
less than zero as z increases, therefore the solution 
in this case is unphysical. 



• Case 3: c < 1 The physical solution is 



H 



1 
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The parameters in this expression should satisfy the 
following condition according to Eq.® at z = 



(11) 



Thereby we obtain JIa and its derivative with respect 
to z in the following forms 



(l-c2)2 
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,e^/flr^ + y/flra + f^m(l " C^) 



3r!,„(l + z) 
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Furthermore, we require that matter dominate over vac- 
uum energy as z grows larger, or else it would spoil the 
success of standard Big Bang cosmology. From Eq. (fT2|) 
we obtain the asymptotic expression of Qa for large z 



nA 



1 



(14) 



Thus if we demand flm dominate over ^Ia at early time, 
we have jz^ <C 1 or c ^ 1/ \/2. In the following we only 
consider the case of c < l/\/2 ^0.7. 



B. L as the particle horizon and the event horizon 

The IR cut-off L is given by the definition of these two 
horizons 



R 



ph 



L = 
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Here we use the Friedmann Eq.Q and recast it into 



H 



E[z) = — ^ 
no 

By Eq.([T|) we have 



L = 
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Inserting Eq.ljTS]) and Eg. p^ into the above equation, 
taking derivative with respect to a on both sides, and 
then using 1 + z = 1/a, we obtain the evolution equation 
of Q,A with respect to z ^ 



n'A = {i+z)-'-nT{ 



where 6 



(18) 



-1 corresponds to L = i?pij and 



-1 



to L = Rch and the initial condition of this differential 
equation is given by setting z=0 in Ea. (|16p 



f2Ao ^1 — ^2^y^l^ — f2 



mO- 



(19) 



We can solve this equation numerically and require that, 
as mentioned above, ^a should become negligible com- 
pared with f2m as z grows. 

In fact, it is only in an eternally accelerating universe 
that the event horizon exists. Thus, when using the event 
horizon as the IR cut-off, we already assume an accelerat- 
ing universe and therefore the existence of some effective 
dark energy. 



C. EoS of the vacuum energy 

By energy conservation we have 

PA + 3H{1 + wa)pa = 0, 
from which we get 

PA 



WA = -1 - 



3ffpA 



(20) 



(21) 



^ It should be noted that here we define the fractional energy den- 
sities with Ho rather than with H, by which these fractions are 
often defined in literature. 



FIG. 1: The evolution of wa{z). L = and fi™ = 0.3. 

In the branch e = +1 (upper), the EoS evolves from zero in 
the past to —1 in the future. It is currently less than —1/3 
and therefore Q.a can serve as dark energy. In the branch 
e = —1 (lower), however, the EoS is always positive. 



FIG. 2: The evolution of WAiz). L — Rch, ~ 0.3 and 
flr^ = 0.12. In both cases the EoS is always less than —1/3 
and the vacuum energy can serve as dark energy. In the 
branch e = +1 (upper), the EoS ends up with —1/3 in the 
future, while in the branch e — —1 (lower) it may cross —1 for 
some values of the parameters. 



where 1 + z = 1/a is used. 

For L as Hubble scale, we insert Eg. p^ and Eq.([T 
into dH]) to get 



-1 + 



(22) 



where c < 1 is required. Some features of the evolution of 
wa can be shown analytically if we rewrite this equation 
as 



-1 + 



1 



e/F^/FTT + F + 1 ' 



(23) 



where F = 



0„(l-c2) 



> 0. For e 



-1, the denominator 



is always greater than 1, leading to wa < forever; for 
e ~ —1, the denominator is always less than 1, leading 
to wa > 0. The evolution of w\ is shown numerically 
in Fig.l, where we fix flmo = 0-3 set different values 
of c. From the figure we find that in the (-) branch the 



EoS is always positive, therefore the vacuum energy can 
not drive the cosmic acceleration. In the (+) branch, 
however, for smaller c ( e.g. c = 0.2,0.4 on the plot), 
w\{z = 0) may become less than —1/3 and the vacuum 
energy may play the role of dark energy. This is different 
from the case in the usual FRW universe, where 
cannot serve as the IR cut-off of the holographic dark 
energy. 

For L as Rph and Rch we insert Eq. pSjl into Eq. lpTjl 
to get 



(24) 



From this equation we can easily see that, it is only when 
6 — —1 that wa < — ^, namely, the vacuum energy serves 
as dark energy if L is R^h rather than Rph- This is the 
same as in the usual FRW universe. Fig[2] shows the 
evolution of the EoS for the case of L as i?oh- 
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D. EoS of the effective dark energy 



In order to explore the possibility of realizing acceler- 
ating expansion in our model, a combined effect of both 
the vacuum energy and the 5D gravity effect should be 
considered. That is, we need to find out the EoS of the 
effective dark energy. Firstly, we rewrite Eg. pB)) as 



E'^ ^ f2„i + fiA + 2f7,.^+2eV^V^m + ^^A + ^r^, (25) 

and then compare this expression with the Friedmann 
equation in usual 4D FRW cosmology consisting of a mat- 
ter component and an effective dark energy 



= fir. 



(26) 



where, by the conservation equation for the effective dark 
energy, we have 



n n(0) /Q r 1 + Wcff(2') , M 
"off = ii^iff exp(3 / — — — - — dz ), (27) 
Jo i + z 

and Wcs is the EoS of the effective dark energy. We find 
that 



(28) 



Taking derivative on both sides of the above equation 
with respect to z, we obtain 

1 1 

1 -I- Wcff = 



(29) 



,^^3a„ + f^:,(i + z) ^ 

Y^m + ^^A + 



with the same constraint as Eg. pop . Here we also require 
the asymptotic behavior of fleS should be dominated over 
by ri„i in the past for the same reason mentioned above. 

For L as the Hubble scale, we show the evolution of 
the effective EoS in FigOand FigH Clearly in the (+) 
branch Wos can become less than —1/3 and end up with 
Weft = — 1, and so is Wtotai due to the effective dark energy 
dominating over matter at late time. For L as the parti- 
cle horizon, FigOand FigHshow that in the (-I-) branch 
the effective dark energy can drive the cosmic accelera- 
tion whereas it cannot in the other branch. Note that, by 
FigEl for small values of c, although the vacuum energy 
may serve as dark energy, it can dominate over matter as 
z grows and consequently spoil the BBN and structure 
formation. As Fig[5]and FigfTUlshow. for L as the future 
event horizon in both branches, as well as Wtotai may 
become less than —1/3. Therefore in both branches ac- 
celeration may occur. Besides, there are two points worth 
particular mentioning: (1) From FigH] and Figl7]we see 
that in the (-) branch, there exist a pole where WcS di- 
verges. This occurs because fles evolves from positive 
to negative (or the opposite) as z decreases, crossing the 
point r^off = at some z*, which implies the breakdown 



-0.2 - 
-0.3 . 




FIG. 3; The evolution of Wcff. L = H~ , e = +1 and = 
0.3. The effective EoS may become less than —1/3 in the 
near past and end up with —1 in the future, therefore an 
acceleration may occur in this case. 



of the effective description, rather than any pathology of 
the model. This can also be confirmed from the plot of 
Wtotai where the EoS is well behaved. (2) The EoS in 
FiglQ] to FiglTT] exhibits phantom behavior. This is be- 
cause the effective dark energy posses phantom behavior 
and it dominates over matter at late time. Therefore fiog 
increases with time until the Big Rip singularity. 



III. PARAMETER FITTING WITH TYPE lA 
SUPERNOVA OBSERVATION AND BARYON 
ACOUSTIC OSCILLATIONS 

In this section, we confront our model with observa- 
tional data and constrain the parameters. We use the 
SNe data compiled by Davis et al.[l3|, which consists 
of 192 SNe classified as SNIa with redshift up to z = 
1.755. This dataset is a combination of several subsets 
which are 45 nearby SNelg], 60 SNe from ESSENCEfTl, 
57 SNe from the Supernova Legacy Survey (SNLS)(2(J| 
and 30 high redshift SNe from Hubble Space Telescope 
(HST)[2l|, of which the data from SNLS and the nearby 
SNe were refitted in What the supernova observa- 
tions provide is the distance modulus fiobj ■ This quantity 
can be calculated from the model by 



Hth{z; {0k}) = 5lgDL{z; {0k}) +M, 



(30) 



where {Ok} represents the parameters of the model: 
{rimOyc} for L = H~^, and {nmo,^r,,c} for L = 
i?ph or i?ch- is a nuisance parameter consisting of 
the Hubble constant Hq and the absolute magnitude M 



(31) 



(here c denotes the speed of light). Dj^ is the dimen- 
sionless luminosity distance free of Hubble constant. In 
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FIG. 4: The evolution of tUcff (upper) and itJtotai (lower). L = 
H-\e=-la.ndQm = 0.3. There is a pole at corresponding 
redshift z* where Weff becomes divergent. This only means 
that the effective description breaks down around z* , rather 
than that some pathology exists in the model. This can be 
illustrated by the plot of TOtotai, which is well behaved. 
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FIG. 6: The evolution of Wcff- L = Rph, e = +1, ~ 0.3 
and fire = 0.12. For c = 0.2 or smaller, lUeff grows larger 
and becomes positive (upper) at high redshift region. This 
implies that the vacuum energy dominates over matter as z 
grows larger (lower), which would violate the success of the 
standard Big Bang cosmology and therefore is not realistic. 
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FIG. 5: The evolution of WcS- L = Rph, e = +1, fim = 0.3 
and fire ~ 0.12. uicff (O) < —1/3 and an acceleration may 
occur. Wcfi — > —1 in the future and the effective dark energy 
ends up as a cosmological constant. 



FIG. 7: The evolution of Wctt- L = Rph- fim = 0.3, fir^ = 
0.12 and e = —1. As the case shown in FIG|4j there exists 
a pole in the future evolution of the EoS, which implies the 
breakdown of the effective description. 
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— c = 0.1 
--c = 0.2 



FIG. 9: The evolution of w^s. L = i?eh, e = +1, f^m = 0.3 
and n,. = 0.06. 



a spatially flat universe it is defined by 



D, 



(1 + ^) 



dz' 





(32) 



The best fits are obtained by minimizing the quantity 
X {{Ok},M) ^ 2.^ , (33) 



where <Ti are the observational uncertainties. We actually 
deal with the quantity with the nuisance parameter 
marginalized over 

= -21n J e-^^l'^dM. (34) 

This is equivalent to minimize x^ with respect to M. [2^ , 
up to a negligible constant. One can easily show that x^ 
can be expanded in M. around = by 

x2({0fc}, M) ^x\M^ 0, {dk}) - 2BM + CM^ (35) 

where 

192 

B = ^ -2 ' (36) 
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where zi = 0.35 is the typical redshift of the LRGs. The 
measured value of A is 0.469±0.017^]. Correspondingly 
the quantity is 



2 _ (A -0.469)2 
Xbao - 



(40) 



0.0172 ■ 

And in the following we perforin the joint analysis using 
X ^Xsn + Xbao- 
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FIG. 12: Contour plot within 3a from a joint analysis for 
L = H'^ and e = +1. The best fits are n^o = 0.25 ± 0.02 
and c = + 0.14. 



FIG. 11: The evolution of Wcb- L = i?ch, e = —1, f^m = 0.3 
and f2rc = 0.06. 



192 

(37) 

i * 

Obviously becomes minimized ii Jv[ — B /C . There- 
fore, in practice we use 

xlNim) ^x'iM^ 0, {e^}) ~ ^ (38) 

as an alternative to for the sake of efficiency in prac- 
tical calculation without losing accuracy. 

Another observational data we resort to as a comple- 
ment to SNe data is from the observations of baryon 
acoustic oscillation peak(BAO). The acoustic oscillations 
in the relativistic plasma at the recombination epoch may 
imprint on the power spectrum of the non-relativistic 
matter of the late universe. And this acoustic signature 
in the large scale clustering of galaxies was detected by 
Eisenstein et al.[23'| using a large sepctroscopic sample of 
the luminous red galaxies (LRGs) from the Sloan Digital 
Sky Survey (SDSS)[23]. We use the model-independent 
parameter A as given in[23| 




c 



FIG. 13: Marginalized likelihood functions for ilmo and c. 



For L = H^^, there are two parameters (r2„io,c) to 
be fitted. And il,.^ is determined by Eg. pTj) . The result 
is shown in Fig[T21 The best fits fimo = 0.25, c = 
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FIG. 14: Contour plot within 3a Irom a joint analysis for 
L = i?ph and e = +1, with c marginalized. The best fits are 
= 0.14 ± 0.01 and Qmo = 0.25 ± 0.02. The dotted line 
denotes the situation of Qao = 0, above which Qao < and 
there the region is forbidden in our model. 




FIG. 15: Marginalized likelihood functions for QmO and Qr^- 



and flr^ — 0.14 indicate that the observations prefer a 
pure DGP model without the holographic vacuum en- 
ergy. This best fits are also consistent with those ob- 
tained in liter atureflBl- Fig [13] shows the marginalized 
likelihood function for the two parameters, in which the 
curve for flmo is near-Gaussian whereas it is highly asym- 
metric for c due to the theoretical cutoff of this parame- 
ter. This leads to flmo = 0.25 ± 0.02 and c < 0.14 within 
68.3% confidence level. 

For L = i?ph, the results are given by FigHJ] and 
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FIG. 16: Contour plot within 3a from a joint analysis for 
L = i?ch and e = -|-1, with c marginalized. The best fits are 
fire =0-1- 0.04 and ^ImO = 0.27 ± 0.02. Corresponding to 
the best fits, c = 0.76. Again, the dotted line denotes the 
boundary between S7ao > O(below) and JIao < O(above). 




FIG. 17: Marginalized likelihood functions for flmo and flr^- 



Fig llSl where we have marginalized the parameter c. In 
the contour plot, the dotted line represents Q\o = in 
Eg. p^ . Below this Hne flAo > 0. In the region above 
^^Ao < and therefore it is the unphysical region for the 
parameters of the model. In the pure DGP model, the 
counterpart of Eq. is 

f^fc = l-2Va~-f^mO, (41) 

where Vlk denotes the spatial curvature. If we identify 
57ao in Eq. (fTO|) with ilk in Eq. (HT1) , we find that our model 
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FIG. 18; Contour plot for Ax^ = 2.30, 6.17, 9.21 from a joint 
analysis for L = Rch and e = — 1, with fimo marginalized. 
The best fits are flr^ =0, c = 0.77 denoted by a star on the 
plot. Corresponding to the best fits, Qmo = 0.27. 



TABLE I; Best fits 

+1 0.25 0.14 -0.8000 
i?ph +1 0.25 0.14 N/A -0.7015 
Roh +1 0.27 0.76 -1.0828 
J?ch -1 0.27 0.77 -1.0731 



in this case is equivalent to the pure DGP model confined 
to a non-closed universe. The best fits are ^Ir^ = 0.14 ± 
0.01 and = 0.25±0.02, indicating a pure DGP model 
in a flat universe without vacuum energy. 

For L = Reh, we have to consider the two branches. 
In the branch, the results are shown in Figlin] and 
FigflTl with the parameter c marginalized. The best fits 
arc ilr^ = + 0.04 and il„io = 0.27 ± 0.02, correspond- 
ingly c = 0.76, indicating a pure holographic dark energy 
model with the negligible effect of higher dimensional 
gravity. In the (-) branch, as we can see from FigfTSl 
(with marginalized), the two outmost contours are 
not closed within a large region of the parameter space, 
indicating that current observations cannot impose tight 
constraint on the parameters in this case. Therefore the 
contours just represent the difference with respect to the 
minimum of x^, = 2.30, 6.17, 9.21 respectively, with- 
out an exact statistical meaning. Thus we do not present 



the likelihood plot as before. Despite of this, we can still 
get the best fits as ilr^ = and c = 0.77 with ilmo = 0.27 
correspondingly. This also indicates a pure holographic 
dark energy model. 

IV. CONCLUSION 

In this paper, we considered the evolution of the vac- 
uum energy in the universe described by the DGP model. 
By numerically studying the EoS of the vacuum energy 
and the effective EoS of the combined effect of both vac- 
uum energy and brane effect, we found that choosing the 
IR cut-off as the event horizon, the vacuum energy can 
drive the cosmic acceleration in both branches. In ad- 
dition, the choice of the Hubble scale as the cut-off can 
also lead to the vacuum energy playing the role of dark 
energy. This is different from the case in ordinary 4D 
gravity, where wa < —1/3 only when the event horizon 
is chosen as the IR cut-off. Further investigation shows 
that when L = i?ch) the EoS may cross —1 and the vac- 
uum energy would end up with a phantom phase, there- 
fore the Big Rip singularity is inevitable, in contrast to 
the models such as LDGP[12j] and SDGPJJ,] where only 
the effective EoS posses the crossing behavior and the 
total EoS is always larger than —1. 

Through a joint analysis of SNe data and BAO data, 
the results of parameter fitting show that the DGP model 
with holographic vacuum energy can be consistent with 
the joint data constraints within 68.3% confidence level. 
For IR cut-off L as the Hubble scale and the particle 
horizon in the (-I-) branch, the best fits indicate that 
the observational data prefer a pure DGP model with 
negligible vacuum energy. For L as the event horizon in 
both branches, on the other hand, the best fits show a 
preference to the pure holographic dark energy model. 
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We consider the evolution of the vacuum energy in the DGP model according to the holographic 
principle under the assumption that the relation linking the IR and UV cut-offs still holds in this 
scenario. The model is studied when the IR cut-off is chosen to be the Hubble scale , the 
particle horizon _Rph and the future event horizon Rc\i, respectively. And the two branches of the 
DGP model are also taken into account. Through numerical analysis, we find that in the cases of 
in the (-f) branch and i?eh in both branches, the vacuum energy can play the role of dark 
energy. Moreover, when considering the combination of the vacuum energy and the 5D gravity 
effect in both branches, the equation of state of the effective dark energy may cross —1, which may 
lead to the Big Rip singularity. Besides, we constrain the model with the Type la supernovae and 
baryon oscillation data and find that our model is consistent with current data within la, and that 
the observations prefer either a pure holographic dark energy or a pure DGP model. 
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I. INTRODUCTION 

Recent SNe la and WMAP observations P, Q have 
indicated that our universe is currently undergoing an 
accelerating expansion, which confront the fundamen- 
tal theories with great challenges and also make the re- 
searches on this problem a major endeavor in modern 
astrophysics and cosmology. The origin of the cosmic ac- 
celeration is still a mystery and is referred to as the dark 
energy problem. Various models have been proposed to 
solve this problem. They generally fall into the following 
two ways. One is to add an exotic energy component 
with negative pressure, that is, the dark energy, to the 
total energy budget of the universe. Among others, the 
most competitive candidate of dark energy sofar is the 
cosmological constant due to both its theoretical simple- 
ness and its great success in fitting with observational 
data, although it suffers from the cosmological constant 
problem^ . Such problem is expected to be solved or al- 
leviated in the models of dynamic dark energy (see for 
a more detailed review), which generally contains a scalar 
field evolving in time and driving the acceleration, just 
like the scalar field introduced for the inflation stage at 
early universe. In fact, the cosmological constant prob- 
lem is essentially a problem of quantum gravity. In quan- 
tum field theory, where the effect of gravity is neglected, 
the vacuum energy is determined by the UV cut-off k^, 
that is PA cx A:^. No matter how we choose the UV cut- 
off, be it the Planck scale lO^^GeV or the electroweak 
scale TeV, the value predicted by theory is far greater 
that that observed IQ-^^GeV*. bmce we are concerning 
problems at the cosmological scale, however, we should 



have take into account the effect of gravity. It is expected 
that the value of the cosmological constant or the vacuum 
energy would be predicted correctly from a complete the- 
ory of quantum gravity, which is still being explored. But 
at present the holographic principle, which is believed to 
be an important feature of quantum gravity, may shed 
some light on solving this problem. Follow the line of 
the holographic principle, the holographic dark energy 
model [7] is a promising candidate for solving the dark 
energy problem. In this model, the vacuum energy is 
no longer a time- independent constant, but evolves with 
time according to the holographic principle. The vacuum 
energy is related to the length measure on cosmological 
scale 
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where pA is directly related to the UV cut-off, L is the 
IR cut-off, Mp denotes the Planck mass and c is a nu- 
merical factor by convention, which is the parameter of 
the model. This relation of the entanglement of UV/IR 
was first proposed in , where L was first chosen as the 
Hubble scale. Then Hsu[§| pointed out that this would 
not lead to the desired equation of state. Finally, LiQ 
proposed the holographic dark energy model where L is 
the event horizon. And this model fits very well with 
current observations [iGj . 

The other way to solve the dark energy problem is 
to modify the theory of gravity at large scale, without 
resorting to any new energy component. For example, 
the /(i?) theorylsll modifies the standard Einstein- Hilbert 
action to introduce an effective dark energy component 
in the Einstein frame. Here we focus on the DGP model 
, which describes our universe as a 4Z? brane embedded 
in a bD Minkowski bulk and explains the origin of the 
dark energy as the gravity on the brane leaking into the 
bulk at large scale. The model is described by the action 
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M2 



(2) 



where gab is the bulk metric and /i^j^ is the induced metric 
on the brane. The first term contains the 51? Ricci scalar 
whereas the second term contains the AD Ricci scalar 
on the brane, which is an extra term due to quantum 
effects, in contrast to the Randall-Sundrum scenario llj . 
The third term represents matter localized on the brane. 
And Sgh is the Gibbsons-Hawking boundary term. 

In this paper, we assume that the relation Eq.([T]) 
still holds in the DGP model, and consider the evo- 
lution of the vacuum energy on the brane (or the 
brane tension) according to the holographic princi- 
ple. Note that there are other models also gener- 
alizing the standard DGP model by adding a cos- 
mological constant (LDGP)[l3|, a Quiessence perfect 
fluid(QDGP)[l3], a scalar field(SDGP)[l|, or the Chap- 
lygin gas(CDGP)[ll. Although the holographic dark en- 
ergy model is well consistent with observational data, it 
should be noted that the core of the holographic principle 
is that it relates the UV and IR cut-offs of a local quan- 
tum field system, which reflects some feature of quantum 
gravity. Thus even if the holographic vacuum energy had 
not played the role of dark energy, it would still be of sig- 
nificance to study on this problem. 

The paper is organized as follows. In section [Til we 
present the model under three cases of choosing the IR 
cut-off for the vacuum energy as respectively the Hubble 
scale H~^, the particle horizon i?ph and the future event 
horizon i?oh- We study the equation of state (EoS) of the 
vacuum energy and the effective EoS due to the combined 
effect of both the holographic vacuum energy and the 5D 
gravity effect. In section Hill we use recent observational 
data to constrain the model and fit its parameters. We 
conclude this paper in the final section. 



(early times), the 4D general relativity is recovered; for 
> (late times), the 5D effect begins to be sig- 
nificant, e = ± represents two branches of the model of 
which the (-I-) branch is the self-accelerating solution in 
which the universe may enter into an accelerating phase 
in late time by virtue of pure 5D effect of gravity, while 
the (— ) branch cannot undergo an acceleration without 
additional dark energy component. Here the vacuum en- 
ergy is added in p = pm + pa- Obviously the dynamics 
is typically different from the standard FRW cosmology. 
In addition, we have the usual equation of conservation 



3H{l + w)p = 0, 



(5) 



where w is the equation of state (EoS). Here we assume 
that there is no interaction between matter and vacuum 
energy. Therefore both components obey the equation of 
conservation respectively and, in particular, for matter 
we have pm = Pm^o,~^ ■, the same as that in usual FRW 
cosmology. We note that here the EoS of the vacuum 
energy evolves with time due to the holographic principle 
as shown later, as opposed to ACDM where wa = — 1. 

When applying the holographic principle in cosmology, 
a crucial problem is how to choose the IR cut-off L. In 
usual FRW sosmology, it is shownQ that only in the case 
of choosing L as the future event horizon can vacuum en- 
ergy play the role of dark energy. In the DGP framework, 
we consider L as the Hubble scale , the particle hori- 
zon i?ph and the future event horizon i?ch respectively. 
And then we investigate the evolution of the EoS of the 
vacuum energy. 

Intuitionally, we may guess that for one thing, in the (-) 
branch where there is no self-acceleration, adding a com- 
ponent whose EoS > —1/3 will not lead to an accelera- 
tion whereas adding a dark energy component may cause 
the universe in this branch to accelerate; for the other 
thing, in the (+) branch, no matter whether the holo- 
graphic vacuum energy we add can in itself play the role 
of dark energy, the combined effect may lead to an accel- 
eration due to the contribution from the self-acceleration 
of this branch. 



II. THE MODEL 



L as H- 



We assume a flat, homogeneous and isotropic brane in 
accordance with the result of the WMAP observation j^]. 
Following 6] , the Friedmann equation is 
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where H =z a/a la the Hubble parameter and Tc = 
Mp/2M^^^ is the distance scale reflecting the competi- 
tion between 4D and 5D effects of gravity. For ^ Vc 



By Eq.([T]), the vacuum energy is p\ = ic^MpH"^. 
For convenience we insert this result into the Friedmann 
Eq.dll). After defining 17„ = = ^™o(l + zf, 



and r^A = 



„2 



where Hq is the Hubble parameter at redshift z — Q, the 
above equation can be transformed into 



nr. 



0, 



(6) 



Ho " Hq 

Then we can solve the above equation to get E = H/Hq 
• Case 1: c = 1 

H 1 



3M2 



(7) 
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where e = +1 corresponds to the contracting sohi- 
tion, and e = — 1 represents the expandmg solution 



H 



3M: 



(8) 



the conservation equation of oc a ^ there- 

fore Pa oc oc (X a^^, namely, the vacuum 



In this case, the solution is easy to find by virtue of 

m ■ Pm ^ ^ 

energy decreases faster than pm and it cannot be 
dominant in late time. Thus we do not consider 
this case as of physical interest. 

When c 7^ 1 we write the general solution 



1-C2 



(9) 



• Case 2: c > 1 Since flr^ is constant, the part under 
the square root on the RHS of Eq. ^ may become 
less than zero as z increases, therefore the solution 
in this case is unphysical. 



• Case 3: c < 1 The physical solution is 



H 



1 



(10) 



The parameters in this expression should satisfy the 
following condition according to Eq.® at z = 



(11) 



Thereby we obtain JIa and its derivative with respect 
to z in the following forms 



(l-c2)2 



(12) 
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,e^/flr^ + y/flra + f^m(l " C^) 



3r!,„(l + z) 



1-C2 
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(13) 



Furthermore, we require that matter dominate over vac- 
uum energy as z grows larger, or else it would spoil the 
success of standard Big Bang cosmology. From Eq. (fT2|) 
we obtain the asymptotic expression of Qa for large z 



nA 



1 



(14) 



Thus if we demand flm dominate over ^Ia at early time, 
we have jz^ <C 1 or c ^ 1/ \/2. In the following we only 
consider the case of c < l/\/2 ^0.7. 



B. L as the particle horizon and the event horizon 

The IR cut-off L is given by the definition of these two 
horizons 



R 



ph 



L = 



Rch - a{t) /^°° ^ _ a XT ^ 



(15) 



Here we use the Friedmann Eq.Q and recast it into 



H 



E[z) = — ^ 
no 

By Eq.([T|) we have 



L = 



VLa -|- f^r 



/3c2Af2 

PA 
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(17) 



Inserting Eq.ljTS]) and Eg. p^ into the above equation, 
taking derivative with respect to a on both sides, and 
then using 1 + z = 1/a, we obtain the evolution equation 
of Q,A with respect to z ^ 



n'A = {i+z)-'-nT{ 



where 6 



(18) 



-1 corresponds to L = i?pij and 



-1 



to L = Rch and the initial condition of this differential 
equation is given by setting z=0 in Ea. (|16p 



f2Ao ^1 — ^2^y^l^ — f2 



mO- 



(19) 



We can solve this equation numerically and require that, 
as mentioned above, ^a should become negligible com- 
pared with f2m as z grows. 

In fact, it is only in an eternally accelerating universe 
that the event horizon exists. Thus, when using the event 
horizon as the IR cut-off, we already assume an accelerat- 
ing universe and therefore the existence of some effective 
dark energy. 



C. EoS of the vacuum energy 

By energy conservation we have 

PA + 3H{1 + wa)pa = 0, 
from which we get 

PA 



WA = -1 - 



3ffpA 



(20) 



(21) 



^ It should be noted that here we define the fractional energy den- 
sities with Ho rather than with H, by which these fractions are 
often defined in literature. 



FIG. 1: The evolution of wa{z). L = H'^ and = 0.3. 
In the branch e = +l(top), the EoS evolves from zero in 
the past to —1 in the future. It is currently less than —1/3 
and therefore JIa can serve as dark energy. In the branch 
e = —1 (bottom), however, the EoS is always positive. 



FIG. 2: The evolution of wa{z). L — Rch, ~ 0.3 and 
flr^ = 0.12. In both cases the EoS is always less than —1/3 
and the vacuum energy can serve as dark energy. In the 
branch e = +l(top), the EoS ends up with —1/3 in the future, 
while in the branch e = —1 (bottom) it may cross —1 for some 
values of the parameters. 



where 1 + z = 1/a is used. 

For L as Hubble scale, we insert Eg. p^ and Eq.([T 
into dH]) to get 



-1 + 



(22) 



where c < 1 is required. Some features of the evolution of 
wa can be shown analytically if we rewrite this equation 
as 



WA 



-1 + 



1 



e/F^/FTT + F + 1 ' 



(23) 



where F = 



0„(l-c2) 



> 0. For e 



-1, the denominator 



is always greater than 1, leading to wa < forever; for 
e ~ —1, the denominator is always less than 1, leading 
to Wa > 0. The evolution of wa is shown numerically 
in Fig.l, where we fix flmo = 0-3 set different values 
of c. From the figure we find that in the (-) branch the 



EoS is always positive, therefore the vacuum energy can 
not drive the cosmic acceleration. In the (+) branch, 
however, for smaller c ( e.g. c = 0.2,0.4 on the plot), 
wa{z = 0) may become less than —1/3 and the vacuum 
energy may play the role of dark energy. This is different 
from the case in the usual FRW universe, where H^^ 
cannot serve as the IR cut-off of the holographic dark 
energy. 

For L as i?ph and Roh we insert Eq. pSjl into Eq. lpTjl 
to get 



WA = -r^ 



3c y^riffi + Ha + 



(24) 



From this equation we can easily see that, it is only when 
6 — —1 that Wa < — ^, namely, the vacuum energy serves 
as dark energy if L is R^h rather than Rph- This is the 
same as in the usual FRW universe. Fig[2] shows the 
evolution of the EoS for the case of L as i?oh- 
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D. EoS of the effective dark energy 



In order to explore the possibility of realizing acceler- 
ating expansion in our model, a combined effect of both 
the vacuum energy and the 5D gravity effect should be 
considered. That is, we need to find out the EoS of the 
effective dark energy. Firstly, we rewrite Eg. pB)) as 



E'^ ^ f2„i + fiA + 2f7,.^+2eV^V^m + ^^A + ^r^, (25) 

and then compare this expression with the Friedmann 
equation in usual 4D FRW cosmology consisting of a mat- 
ter component and an effective dark energy 



= fir. 



(26) 



where, by the conservation equation for the effective dark 
energy, we have 



n n(0) /Q r 1 + Wcff(2') , M 
"off = ii^iff exp(3 / — — — - — dz ), (27) 
Jo i + z 

and Wcs is the EoS of the effective dark energy. We find 
that 



(28) 



Taking derivative on both sides of the above equation 
with respect to z, we obtain 

1 1 

1 -I- Wcff = 



(29) 



,^^3a„ + f^:,(i + z) ^ 

Y^m + ^^A + 



with the same constraint as Eg. pop . Here we also require 
the asymptotic behavior of fleS should be dominated over 
by ri„i in the past for the same reason mentioned above. 

For L as the Hubble scale, we show the evolution of 
the effective EoS in FigOand FigH Clearly in the (+) 
branch Wos can become less than —1/3 and end up with 
Weft = — 1, and so is Wtotai due to the effective dark energy 
dominating over matter at late time. For L as the parti- 
cle horizon, FigOand FigHshow that in the (-I-) branch 
the effective dark energy can drive the cosmic accelera- 
tion whereas it cannot in the other branch. Note that, by 
FigEl for small values of c, although the vacuum energy 
may serve as dark energy, it can dominate over matter as 
z grows and consequently spoil the BBN and structure 
formation. As Fig[5]and FigfTUlshow. for L as the future 
event horizon in both branches, as well as Wtotai may 
become less than —1/3. Therefore in both branches ac- 
celeration may occur. Besides, there are two points worth 
particular mentioning: (1) From FigH] and Figl7]we see 
that in the (-) branch, there exist a pole where WcS di- 
verges. This occurs because fles evolves from positive 
to negative (or the opposite) as z decreases, crossing the 
point r^off = at some z*, which implies the breakdown 



-0.2 - 
-0.3 . 




FIG. 3; The evolution of Wcff. L = H~ , e = +1 and = 
0.3. The effective EoS may become less than —1/3 in the 
near past and end up with —1 in the future, therefore an 
acceleration may occur in this case. 



of the effective description, rather than any pathology of 
the model. This can also be confirmed from the plot of 
Wtotai where the EoS is well behaved. (2) The EoS in 
FiglQ] to FiglTT] exhibits phantom behavior. This is be- 
cause the effective dark energy posses phantom behavior 
and it dominates over matter at late time. Therefore fiog 
increases with time until the Big Rip singularity. 



III. PARAMETER FITTING WITH TYPE lA 
SUPERNOVA OBSERVATION AND BARYON 
ACOUSTIC OSCILLATIONS 

In this section, we confront our model with observa- 
tional data and constrain the parameters. We use the 
SNe data compiled by Davis et al.[l3|, which consists 
of 192 SNe classified as SNIa with redshift up to z = 
1.755. This dataset is a combination of several subsets 
which are 45 nearby SNelg], 60 SNe from ESSENCEfTl, 
57 SNe from the Supernova Legacy Survey (SNLS)(2(J| 
and 30 high redshift SNe from Hubble Space Telescope 
(HST)[2l|, of which the data from SNLS and the nearby 
SNe were refitted in What the supernova observa- 
tions provide is the distance modulus fiobj ■ This quantity 
can be calculated from the model by 



Hth{z; {0k}) = 5lgDL{z; {0k}) +M, 



(30) 



where {Ok} represents the parameters of the model: 
{rimOyc} for L = H~^, and {nmo,^r,,c} for L = 
i?ph or i?ch- is a nuisance parameter consisting of 
the Hubble constant Hq and the absolute magnitude M 



(31) 



(here c denotes the speed of light). Dj^ is the dimen- 
sionless luminosity distance free of Hubble constant. In 
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FIG. 4: The evolution of i(;cff(top) and ?Utotai (bottom). L = 
e = -1 and Qm = 0.3. There is a pole at corresponding 
redshift z* where Weff becomes divergent. This only means 
that the effective description breaks down around z* , rather 
than that some pathology exists in the model. This can be 
illustrated by the plot of lutotai, which is well behaved. 
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FIG. 6: The evolution of Wctf. L = Rph, e = +1, fim ~ 0.3 
and r^jc ~ 0.12. For c — 0.2 or smaller, WcS grows larger and 
becomes positive (top) at high redshift region. This implies 
that the vacuum energy dominates over matter as z grows 
larger (bottom), which would violate the success of the stan- 
dard Big Bang cosmology and therefore is not realistic. 
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FIG. 5: The evolution of WcS- L = Rph, e = +1, fim = 0.3 
and fire ~ 0.12. uicff (O) < —1/3 and an acceleration may 
occur. Wcfi — > —1 in the future and the effective dark energy 
ends up as a cosmological constant. 



FIG. 7: The evolution of Wctt- L = Rph- = 0.3, D,r^ = 
0.12 and e = —1. As the case shown in FIG|4j there exists 
a pole in the future evolution of the EoS, which implies the 
breakdown of the effective description. 
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— c = 0.1 
--c = 0.2 



FIG. 9: The evolution of w^s. L = i?eh, e = +1, f^m = 0.3 
and n,. = 0.06. 



a spatially flat universe it is defined by 



D, 



(1 + ^) 



dz' 





(32) 



The best fits are obtained by minimizing the quantity 
X {{Ok},M) ^ 2.^ , (33) 



where <Ti are the observational uncertainties. We actually 
deal with the quantity with the nuisance parameter 
marginalized over 

= -21n J e-^^l'^dM. (34) 

This is equivalent to minimize x^ with respect to M. [2^ , 
up to a negligible constant. One can easily show that x^ 
can be expanded in M. around = by 

x2({0fc}, M) ^x\M^ 0, {dk}) - 2BM + CM^ (35) 

where 

192 

B = ^ -2 ' (36) 



8 



20 






10 













-10 






-20 






-30 












— c = 0.1 

- - c = 0.2 



-0.5 0.5 



-c = 0.1 
-C = 0.2 



where zi = 0.35 is the typical redshift of the LRGs. The 
measured value of A is 0.469±0.017^]. Correspondingly 
the quantity is 



2 _ (A -0.469)2 
Xbao - 



(40) 



0.0172 ■ 

And in the following we perforin the joint analysis using 
X ^Xsn + Xbao- 
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FIG. 12: Contour plot within 3a from a joint analysis for 
L = H'^ and e = +1. The best fits are n^o = 0.25 ± 0.02 
and c = + 0.14. 



FIG. 11: The evolution of Wcb- L = i?ch, e = —1, f^m = 0.3 
and f2rc = 0.06. 



192 

(37) 

i * 

Obviously becomes minimized ii Jv[ — B /C . There- 
fore, in practice we use 

xlNim) ^x'iM^ 0, {e^}) ~ ^ (38) 

as an alternative to for the sake of efficiency in prac- 
tical calculation without losing accuracy. 

Another observational data we resort to as a comple- 
ment to SNe data is from the observations of baryon 
acoustic oscillation peak(BAO). The acoustic oscillations 
in the relativistic plasma at the recombination epoch may 
imprint on the power spectrum of the non-relativistic 
matter of the late universe. And this acoustic signature 
in the large scale clustering of galaxies was detected by 
Eisenstein et al.[23'| using a large sepctroscopic sample of 
the luminous red galaxies (LRGs) from the Sloan Digital 
Sky Survey (SDSS)[23]. We use the model-independent 
parameter A as given in[23| 




c 



FIG. 13: Marginalized likelihood functions for ilmo and c. 



For L = H^^, there are two parameters (r2„io,c) to 
be fitted. And il,.^ is determined by Eg. pTj) . The result 
is shown in Fig[T21 The best fits fimo = 0.25, c = 
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FIG. 14: Contour plot within 3a Irom a joint analysis for 
L = i?ph and e = +1, with c marginalized. The best fits are 
= 0.14 ± 0.01 and Qmo = 0.25 ± 0.02. The dotted line 
denotes the situation of Qao = 0, above which Qao < and 
there the region is forbidden in our model. 




FIG. 15: Marginalized likelihood functions for QmO and Qr^- 



and flr^ — 0.14 indicate that the observations prefer a 
pure DGP model without the holographic vacuum en- 
ergy. This best fits are also consistent with those ob- 
tained in liter atureflBl- Fig [13] shows the marginalized 
likelihood function for the two parameters, in which the 
curve for flmo is near-Gaussian whereas it is highly asym- 
metric for c due to the theoretical cutoff of this parame- 
ter. This leads to flmo = 0.25 ± 0.02 and c < 0.14 within 
68.3% confidence level. 

For L = i?ph, the results are given by FigHJ] and 



0.2'.. 
0.18 - 
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0.14 . 
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FIG. 16: Contour plot within 3a from a joint analysis for 
L = i?ch and e = -|-1, with c marginalized. The best fits are 
fire =0-1- 0.04 and ^ImO = 0.27 ± 0.02. Corresponding to 
the best fits, c = 0.76. Again, the dotted line denotes the 
boundary between S7ao > O(below) and JIao < O(above). 




FIG. 17: Marginalized likelihood functions for flmo and flr^- 



Fig llSl where we have marginalized the parameter c. In 
the contour plot, the dotted line represents Q\o = in 
Eg. p^ . Below this Hne flAo > 0. In the region above 
^^Ao < and therefore it is the unphysical region for the 
parameters of the model. In the pure DGP model, the 
counterpart of Eq. is 

f^fc = l-2Va~-f^mO, (41) 

where Vlk denotes the spatial curvature. If we identify 
57ao in Eq. (fTO|) with ilk in Eq. (HT1) , we find that our model 
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FIG. 18; Contour plot for Ax^ = 2.30, 6.17, 9.21 from a joint 
analysis for L = Rch and e = — 1, with fimo marginalized. 
The best fits are flr^ =0, c = 0.77 denoted by a star on the 
plot. Corresponding to the best fits, Qmo = 0.27. 



TABLE I; Best fits 

+1 0.25 0.14 -0.8000 
i?ph +1 0.25 0.14 N/A -0.7015 
Roh +1 0.27 0.76 -1.0828 
J?ch -1 0.27 0.77 -1.0731 



in this case is equivalent to the pure DGP model confined 
to a non-closed universe. The best fits are ^Ir^ = 0.14 ± 
0.01 and = 0.25±0.02, indicating a pure DGP model 
in a flat universe without vacuum energy. 

For L = Reh, we have to consider the two branches. 
In the branch, the results are shown in Figlin] and 
FigflTl with the parameter c marginalized. The best fits 
arc ilr^ = + 0.04 and il„io = 0.27 ± 0.02, correspond- 
ingly c = 0.76, indicating a pure holographic dark energy 
model with the negligible effect of higher dimensional 
gravity. In the (-) branch, as we can see from FigfTSl 
(with marginalized), the two outmost contours are 
not closed within a large region of the parameter space, 
indicating that current observations cannot impose tight 
constraint on the parameters in this case. Therefore the 
contours just represent the difference with respect to the 
minimum of x^, = 2.30, 6.17, 9.21 respectively, with- 
out an exact statistical meaning. Thus we do not present 



the likelihood plot as before. Despite of this, we can still 
get the best fits as ilr^ = and c = 0.77 with ilmo = 0.27 
correspondingly. This also indicates a pure holographic 
dark energy model. 

IV. CONCLUSION 

In this paper, we considered the evolution of the vac- 
uum energy in the universe described by the DGP model. 
By numerically studying the EoS of the vacuum energy 
and the effective EoS of the combined effect of both vac- 
uum energy and brane effect, we found that choosing the 
IR cut-off as the event horizon, the vacuum energy can 
drive the cosmic acceleration in both branches. In ad- 
dition, the choice of the Hubble scale as the cut-off can 
also lead to the vacuum energy playing the role of dark 
energy. This is different from the case in ordinary 4D 
gravity, where wa < —1/3 only when the event horizon 
is chosen as the IR cut-off. Further investigation shows 
that when L = i?ch) the EoS may cross —1 and the vac- 
uum energy would end up with a phantom phase, there- 
fore the Big Rip singularity is inevitable, in contrast to 
the models such as LDGP[12j] and SDGPJJ,] where only 
the effective EoS posses the crossing behavior and the 
total EoS is always larger than —1. 

Through a joint analysis of SNe data and BAO data, 
the results of parameter fitting show that the DGP model 
with holographic vacuum energy can be consistent with 
the joint data constraints within 68.3% confidence level. 
For IR cut-off L as the Hubble scale and the particle 
horizon in the (-I-) branch, the best fits indicate that 
the observational data prefer a pure DGP model with 
negligible vacuum energy. For L as the event horizon in 
both branches, on the other hand, the best fits show a 
preference to the pure holographic dark energy model. 
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